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Abstract 

Within the standard model (SM), the time dependent CP asymmetries in B —* 
ipKg, B — > n'K s and B — > cpK s are expected to give the same result i.e. sin 2(3. 
However, recent measurements of the mixing induced CP asymmetries in B — > 4>Kg 
and B — > n'K s modes give results whose central values differ from the SM expectations. 
We explore the effect of new physics in the two Higgs doublet model (THDM), which 
allows tree level flavor changing neutral currents (so called Model III), and the model 
with extra vector-like down quark (VLDQ). We find that the observed mixing induced 
CP asymmetry for B — > <j)Ks can not be accommodated by the THDM, but can be 
explained in the VLDQ model and both models can explain the observed asymmetry 
for B — > n'Ks mode. 



1 Introduction 



A new era in B-physics has just been started with the advent of 5-factories. With the 
accumulation of huge data in the B- system, the standard model (SM) will be subjected to a 
very stringent test. At the same it is also considered that the experiments at B factories are 
also potential sources for probing new physics. The BaBar [1] and Belle [2] measurements of 
the time-dependent asymmetries in the gold plated mode B — > ipKs have provided the first 
evidence of CP violation in the £>-system. The observed world average of sin 2/5 [3], 



agrees well with the SM prediction. This indicates that CP symmetry is significantly violated 
in nature and the Kobayashi-Maskawa (KM) mechanism [4] seems to be the dominant source 
of CP violation, in which the phase 5 KM is the only source of CP violation. However, this 
speculation does not exclude interesting CP violating new physics (NP) effects in other 
B decays. Since the decay B — > ipK$ (b — > ccs) is a tree level process in the SM, the NP 
contributions to its amplitude are naturally suppressed. Moreover, at loop level NP may give 
large contributions to the B°-B° mixing as well as to the loop-induced decay amplitudes. 
The former effects are universal to all -B decay modes and are constrained to be less than 
20% compared to that of the SM contribution [3]. On the other hand, the effects of new 
physics in the decay amplitudes are non- universal, may vary from process to process and 
can show up in the comparison of the CP asymmetries in different decay modes [5]. 

One of the most promising processes for NP searches widely considered in literature [6- 
14] is the decay B — > 4>K S . Various NP scenarios have been presented to explain the data. 
Unfortunately, we do not know at present which is the correct one. Hopefully, careful study 
in future will rule out some of the scenarios, at least as far as the understanding of B physics 
and CP violation is concerned. 

Unlike B — > ipK$, the process B — > <pKs has no tree level amplitude, which makes 
inroads for NP to play an important role in this mode. In the SM the decay b — > sss, 
which contributes to B — > (pKs, is induced at one loop level. Thus, it is natural to expect 
that new physics contribution to this decay mode may be quite significant. According to 
the KM mechanism of CP violation, both CP asymmetries in B — > (j)Kg and B — > ipK$ 
processes should measure the same quantity, namely sin 2/5, with negligible hadronic uncer- 
tainties (upto 0(\ 2 ), A w 0.2) [5, 12]. However, contrary to the SM expectations, the recent 
measurements of CP asymmetries in B — > <pKs by BaBar [15] and Belle [2] collaborations 
have registered significant deviation from the predictions, as 



sin2/Vs = 0.734 ±0.054, 



(1) 



sm(2(3) 4>Ks = -0.191^ ±0.09 BaBar 
sin(2/5)^ Xs = -0.73 ±0.64 ±0.18 Belle 



(2) 



with an average 



sin ( (2/5) fKg) 



-0.39 ±0.41 



(3) 
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The corresponding branching ratio is given (in units of 10 6 ) as 



BR{B - 




= 8.7^ ±0.9 


BaBar 


BR(B - 


+ <f>K°) 


- inn+ L9 +°- 9 

— iu.u_ 17 _ 13 


Belle 


BR(B - 


+ <f>K°) 


= 5.4^ ±0.7 


CLEO [16] 



with an average 

BR(B -> (pK°) ave = 8.67 ± 1.28 . (5) 

One can see that there are large statistical errors associated with these measurements. 
Nevertheless, the data establish a 2.7a deviation from the SM prediction sin(2j3) ( pK s = 
sin(2P) t pK s - Therefore, if the measurement of sin 2/3 in B — > ipKs is considered as the 
first evidence of large CP violation in B system then the difference between sm(2(3) ( p Ks and 
sin(2j3)^K s i s likely to be regarded as a potential hint for the presence of new physics. There 
are several attempts in the literature [6-14] with detail discussion on the possible implications 
of this result. 

The second channel we are interested in is B° — > i]'K S - This is another two-body decay 
mode which is similar to the two mentioned above. Since many alternative schemes have 
been presented in the literature recently to explain the sin(2/3) ( f, Ks deviation, it is therefore 
very important to verify that each of the NP scenarios should successfully explain them all. 
At present it is difficult to say which is the correct description. In order to narrow down the 
same it is highly desirable that one should carefully study them. This will not only help us 
to narrow down the sources of NP but also provide important clues for hadronic B-physics 
in general. 

B° — > 7]'Ks also receives dominant contribution from the b — ^ sss gluonic penguin, and 
therefore it is expected that the time dependent mixing induced CP asymmetry for this 
mode will also give the value sin 2f3 [3] . However, this decay mode also has a tiny CKM as 
well as color suppressed b — > uus tree contributions along with b — > sqq (q = u, d) penguins, 
which induce deviation from the leading result. It has been shown in Ref. [17] that this 
deviation will be below two percent level. Belle [18] and BaBar [19] collaborations have 
recently measured the CP asymmetry for this mode which is given as 

sin(2/3)^ s = 0.71 ± 0.371^6 Belle 

sm(2p) v , Ks = 0.02 ±0.34 ±0.03 BABAR (6) 

with an average 

sin((2/?V*r s ) „ e = 0.33 ±0.25, (7) 

whose central value also deviates significantly from SM expectations. 

In this paper we would like to investigate the new physics effects on the CP asymmetry 
parameters of the decay B° — > <pK s and B° — > rj'K s modes, arising from some simple 
extensions of the SM. The models considered here are the two Higgs doublet models (THDM) 
which allows tree level flavor changing neutral currents, the so called model III and the model 
with extra vector-like down quarks (VLDQ). We show that the observed data for B° — > 4>Ks 
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can be easily accommodated in the VLDQ model whereas it can not be explained in the 
THDM, and both the models can explain the data for B° — > rj'Ks mode. It has already been 
discussed in Ref. [6], that whether these two models can explain the observed CP asymmetry 
in B — > <f)Ks mode, i.e. sm(2(3) ( / ) K s - However in this paper we have explicitly done the 
calculation for both the decay modes and confirm the result of Ref. [6] for the decay mode 
B^<PK S . 

The paper is organized as follows. In section II, we present the basic formulae for CP 
violating parameters, in the presence of new physics. In section III, we discuss CP violation 
effects in B° — > <pKs mode arising from the THDM and VLDQ model. The B° — > rj'Ks 
process is discussed in section IV. Section V contains our conclusion. 



2 CP violation parameters 

Here, we will present the basic formulae of CP asymmetry parameters, in the presence of new 
physics. Due to the contributions from new physics, these parameters deviate substantially 
from their standard model values. Let us consider the B° and B° decay into a CP eigenstate 
fcp ( we consider f CP = <pK s or rj'K s ). Here, we are presenting the formulae for B° — > 4>K S 
mode, but the same results will also hold for B° — > rj'Ks mode. The time dependent CP 
asymmetry for B — > (pK s can be described by [20] 

Atpxs (t) = C 4>Ks cos(AM Bd t) + S^s sm(AM Bd t) (8) 

where we identify 

C i-'Aj 2 _ 2Im(A) 

as the direct and the mixing-induced CP asymmetries. The parameter A corresponds to 

q AjB^Ks) ( . 

X -pA(B^Wi) (10) 

where, q and p are the mixing parameters and represented by the CKM elements in the 
standard model as 

v v tb v* d 1 ) 

Using CKM unitarity, the amplitude for B — ► <J)K S is given as [12, 21] 

A(B -> cjyKs) = X C A CS + \ U A US (12) 

where X q = V g bV* s . The first term which is the dominant one, is real. Thus if one neglects 
the subdominant amplitude i.e. the doubly Cabibbo supressed second term which in general 
expected to be very small, the mixing induced CP asymmetry is given as, S^Ks = s ^ n 2/3, 
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same as the one for B — > ipK$ in the SM. It has beeen shown in Ref. [12], that the correction 
due to the second term is upto 0(\ 2 ) i.e., 

| 1 S 0Xs -sin2/3|<O(A 2 ). (13) 

Adding a mild dynamical assumption to the SU(3) analysis, recently it has been shown in 
Ref. [21] that the upper bound of standard model pollution to the dominant amplitude of 
B — > <f)Ks mode is of the order of 0.25 and for B — > rj'Ks as 0.3. 

New Physics could in principle contribute to both mixing and decay amplitudes. The new 
physics contribution to mixing is universal while it is non universal and process dependent in 
the decay amplitudes. As the NP contributions to mixing phenomena is universal, it will still 
set S^Ks = S<t>K s - Therefore, to explain the observed 2.7a deviation in (S^k s ~ S<t>K s ), nere 
we explore the NP effects only in the decay amplitudes. Thus including the NP contributions, 
we can write the decay amplitude for B — > (j)K process as 



A(B° - <j>K) = Asm + A NP = A SM [1 + r NP e^j (14) 

where r^p = \Anp/Asm\, (Asm and A^p correspond to the SM and NP contributions to 
the B — > (j)K s decay amplitude) and 4>np = A-Yg(A NP / 'Asm), which contains both strong 
and weak phase components. 

The branching ratio for B — > (j)K decay process can be given as 

BR(B -> <f>K) = BR SM (l + r 2 NP + 2r NP cos <p NP ) (15) 



where BR SM represents the corresponding standard model value. 

Now if we write 4>np = Snp + &np, where 5np and 6np are the relative strong and weak 
phases between the new physics contributions to the decay amplitude and the SM part, one 
can then obtain the expressions for the CP asymmetries as 

sm2/3 + 2r NP cos S NP sm(2/3 + 9 NP ) + r 2 vp sm(2/3 + 26 NP ) 

^<t> K ~ 1 — \ — 2 7~^> r a v 1 "/ 

1 + r z NP + 2r NP coso NP cosO NP 

and 

-2r NP sin 5 NP sin 9 NP 

^4>K — — — 2 x a 

1 + rfjp + 2tnp cos onp cos&np 

In Eqs. (16) and (17) there are three unknowns, namely, r^p, 6np and 5np- So if 
somehow we could constrain the value of rj^p considering different new physics models, we 
could vary the 6np and Snp parameters to obtain the required value of S^k- 



3 CP Violation in B —> (j)Ks process 

To study the CP violation effects in B° — > <pKs process, first we present the SM amplitude 
and then we consider the THDM and thereafter the model with extra vector-like down quark, 
in the following subsections. 
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3.1 SM contributions 



In the SM, the decay process B — > <pKs proceeds through the quark level transition b — > sss, 
which is induced by the QCD, electroweak and magnetic penguins. QCD penguins with 
the top quark in the loop contribute predominantly to such process. However, since we are 
looking for NP here we would like to retain all the contributions. The effective Hamiltonian 
describing the decay b — > sss [22] is given as 
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Heff = ( E CPi + C 9 O g ] , (18) 



U'=3 



where 3 , • • • , 6 and 7 , • • • , O w are the standard QCD and EW penguin operators, re- 
spectively, and O g is the gluonic magnetic operator. Within the SM and at scale M w , the 
Wilson coefficients C±(M W ), ■ ■ ■ ,C W (M W ) at next to leading logarithmic order (NLO) and 
C g {M w ) at leading logarithmic order (LO) have been given in Ref [23]. The corresponding 
QCD corrected values at the energy scale \i = m^, can be obtained using the renormalization 
group equation, as described in Ref. [24]. 

To calculate the B meson decay rate, we use the factorization approximation to evaluate 
the hadronic matrix element (Oj) = (K°<f)\Oi\B ) . Since the hadronic matrix elements do 
not appear in the expressions for CP asymmetry parameters, they will not introduce any 
uncertainties in the results. In this approximation the matrix elements are given as (0 3 ) = 
(0 4 ) = AX/3, (0 5 > = X, (O e ) = X/3, (0 7 ) = -X/2, (0 8 > = -X/6 and (0 9 ) = (O 10 ) = 
— 2X/3, where the factorizable hadronic matrix element X is given as X = (X (p^)|s7 jU (l — 
75 )6|5 o (p B ))(0(g,e )|s 7 ^(l - 75 )s|0) = 2F? K (M*)f 4> M 4> e (r p K . For evaluating the matrix 
element of the most relevant operator, i.e., O g , we use the procedure of [25], where it has been 
shown that the operator O g is related to the matrix element of the QCD and electroweak 
penguin operators as 



(Og 



4:71 



(<1 2 ) 



(o 4 ) + (o 6 ) - ^«o 3 ) + (0 5 )) 



(19) 



where g M is the momentum transferred by the gluon to the (s, s) pair. The parameter (q 2 ) 
introduces certain uncertainty into the calculation. In the literature its value is taken in 
the range 1/4 < (q 2 )/ml < 1/2 [26], and we will use (q 2 )/ml = 1/2 [24], in our numerical 
calculations. 

Thus, in the factorization approach the amplitude A = (<pK°\H e ff\B ) of the decay 
B° — > 4>K° takes a form 



A(B° - 0K°) = -% V tb V t 



V2 



is 



a 3 + a 4 + a 5 - -(a 7 + a 9 + ai 0y 



X , (20) 



where X stands for the factorizable hadronic matrix element of which exact form is irrelevant 
for us since it cancels out in the CP asymmetries. The coefficients Oj are given by 

G2i-i = Cll-i + J^~02i f ; a>2i = C^J + j-j-C^-I-L , (21) 
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where N c is the number of colors. The values of the QCD improved effective coefficients 
C\ can be found in [24, 27]. Now substituting the values of Oj for Nc=3, from [27], the 
value of the form factor Ff K (M^) = 0.39 and using the <fi meson decay constant = 0.233 
GeV and tb° = 1.542 x 10~ 12 sec [28], we obtain the branching ratio in SM as 

BR SM (B -> (/)K ) = 10.5 x 10~ 6 (22) 

which lies within the present experimental limits (5). 



3.2 Two Higgs Doublet Model contributions 

We now proceed to calculate the new physics effect in two Higgs doublet Model (THDM), 
which is one of the simplest extensions of the SM [29]. In such models, the tree level flavor 
changing neutral currents (FCNC's) are prevented by imposing one ad hoc discrete symmetry 
to constrain the THDM scalar potential and Yukawa Lagrangian and thus one obtains the 
so called model I and model II [30] . In model I both the up and down type quarks get mass 
from the Yukawa couplings to the same Higgs doublet 4>i an d in Model II the up- and down 
type quarks get their masses from Yukawa couplings to two different scalar doublets 0i and 
02- Here we consider the model III [31] of THDM where no discrete symmetry is imposed 
and both up- and down type quarks may have diagonal and/or off diagonal flavor changing 
couplings with the two Higgs doublets 0i and 2 - 

The Yukawa Lagrangian of the quarks in model III is given in the form [27] 

4" = V^QiJiU^R + vSQi^iDj,R + iYjQi,LhUj,R + iSQi^D^R + h.c (23) 

where <pi (i = 1,2) are the two Higgs doublets of THDM, 0, = ir 2 0*, Q i: L with (i = 1,2,3) 
are the left handed isodoublet quarks U^r (Dj R ) are the right handed isosinglet up (down) 
type quarks. T^j correspond to the diagonal mass matrices of the up and down quarks, 
while the neutral and charged flavor changing couplings are 

Sij — ~ ^neutral ~ S ? ^charged ~ S V CKM, ^neutral ~ V CKMC, , {^) 

where Vqkm is the Cabibbo-Kobayashi-Maskawa mixing matrix [4] . The coupling constants 
Xij are the free parameters of the model to be determined from experimental data. 

Recently Chao et al [32] studied the b — > S7 process and Xiao et al [27] studied the 
charmless nonleptonic decays of B mesons using the model III of THDM where they have 
kept only the couplings \ tt = \\tt\e tet and A&& = | A&&| e 6 as nonzero. From the studies of 
[27, 32], it is known that, the following parameter space for model III 

Xij = for ij ^ tt or bb 

\\ tt \ = 0.3, |A W | = 35, 9 = (0° - 30°), M H+ = (200 ± 100) GeV (25) 
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where 9 = 6 b — 6 t , are allowed by the available data. The advantage of keeping only these 
two couplings nonzero is that the neutral Higgs boson do not contribute at the tree level or 
one loop level. The new contributions therefore come only from the charged Higgs penguin 
loop with heavy internal top quark. 

The new physics will manifest itself by modifying the corresponding Inami-Lim [33] func- 
tions C (x), D (x), E (x) and E' (x) which determine the Wilson coefficients C 3 (M W ), ■ ■ ■ , C W (M W ) 
and Cg(Mw) in SM. The new strong and electroweak penguin diagrams in THDM can be 
obtained from the corresponding penguin diagrams in SM by replacing the internal W ± lines 
by the charged Higgs H ± lines. Following the same procedure as in the SM, it is straight 
forward to calculate the new 7-, Z Q and gluonic penguin diagrams induced by the exchange 
of charged Higgs bosons in Model III. These new Wilson coefficients C^ I± (M W ) % — 3, • • • , 10 
at NLO level and C g at the LO level can now be written as 



cf{M w 
cf{M w 
ct(Mw 

C^(M W 
Ct(M w 

cf{M w 



a s (M w ) NP a em 



24tt 
®s(M w )_ eNP 



(j 

Qn sin 2 8w ° 



NP 



8tt 
a s (M w ) 

24tt 
a s (M w )_ ENP 



E, 



NP 



in 



a, 



6n 



NP 



C? (M w ) = 



67T 

1 



AC^ P + D^ + - 



NP 



-AC, 



NP 



sin 2 9 W ^° 



E 



NP 



(26) 



where the functions Cq P , Dq P , Eq P and E' NP are the new physics contributions to the 
Wilson coefficients arising from the charged Higgs exchange penguin diagrams. These are 
given by 



r<NP 



jjNP 



E 



E 



NP 


'NP 




Xt_ 

16 



Vt 



Vt 



l-Vt 



yt)' 



In y t 



~H(y t )\\t\" 
-^(y t )|A tt | 2 

l -J{y t )\\t\ 2 -K{y t )\\ tt \ bb \e*' 
6 



(27) 
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with 



^ = 72(1 - y) 3 + 12(1 - „)- ^ 
u I69 - 29j/ 2 + 7J/ 3 2j/ - 3y 2 

= 36(1 -„)» + *KT^r k!/ 

% + 5j/ 2 - j/ 3 3j/ 2 , 

*<»> = w=w~ w^w ln s ' (28) 

In the above use has been made of x t = mf/M^ and y t = mf/Mf I+ . 

Since the charged Higgs bosons appeared in Model III have been integrated out at the 
scale M w , the QCD running of Wilson coefficients C^ I± {M W ) down to the scale (j, = 0(m b ) 
using the renormalization group equation can be done in the same way as in the SM. Including 
the new physics contributions the values of the effective Wilson coefficients at the scale Oijn^) 
are explicitly given in Ref. [27]. Using the values for the Wilson coefficients from [27], we 
obtain the B -> <pK° amplitude in THDM as 

A THDM (B -> (f>K ) = A SM (l + 0.28 e i{dNP+&Np) ) . (29) 

Now taking r/yp = 0.28 and varying the weak phase 6>/vp = {— 7r,7r} and strong phase 
5jvp = {0,27r} according to Eq. (16), we find that the value of S^k can n °t be negative as 
shown in Fig-1. Thus the observed value of S^k can n °t be accommodated in the THDM. 



3.3 Contributions from Model with extra vector like down quark 

Now we consider the model with an additional vector like down quark [34]. It is a simple 
model beyond the SM with an enlarged matter sector with an additional vector-like down 
quark D 4 . The most interesting effects in this model concern CP asymmetries in neutral B 
decays into final CP eigenstates. At a more phenomenological level, models with isosinglet 
quarks provide the simplest self consistent framework to study deviations of 3 x 3 unitarity 
of the CKM matrix as well as, allow flavor changing neutral currents at the tree level. The 
presence of an additional down quark implies a 4 x 4 matrix Vi a (i = u, c, t, 4, a — d, s, b, b'), 
diagonalizing the down quark mass matrix. For our purpose, the relevant information for 
the low energy physics is encoded in the extended mixing matrix. The charged currents are 
unchanged except the Vqkm is now the 3x4 upper submatix of V. However, the distinctive 
feature of this model is that FCNC enters neutral current Lagrangian of the left handed 
downquarks : 
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2 cos 9 W 



uu-fuu - d La U aP ^d L(3 - 2 sin 2 9 w J» m (30) 
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0np 90 

18 0° 



Figure 1: 3DPlot of S^k versus the weak phase 6^p and strong phase 5np (in degrees) for 
r NP = 0.28 

with 

Uafi= E V^ = 8 aP -V: a V ip (31) 

i=u,c,t 

where U is the neutral current mixing matrix for the down sector which is given above. As 
V is not unitary, U ^ 1. In particular its non-diagonal elements do not vanish : 

U a p = -V; a V 4/3 ^ fora^/5 (32) 

Since the various U a p are nonvanishing they would signal new physics and the presence 
of FCNC at tree level, this can substantially modify the predictions for CP asymmetries. 
The new element U sb which is relevant to our study is given as 

U sb = V: s V ub + V c * s V cb + V*V tb (33) 

The decay mode B° — ► <pKs receives the new contributions both from color allowed and 
color suppressed Z-mediated FCNC transitions. The new additional operators are given as 

Of- FCNC = [s a Y(l - l 5 )b a ][s^(C v - C s Alb )s p ] 

Ot FCNC = [SffQ. - l,)b a ][s al ^C a v - C' A ^)ap\ (34) 
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where C v and C% are the vector and axial vector Zss couplings. Using Fierz transformation 
and the identity (C s v - C% lb ) = [{C s v + C%)(1 - 75) + (C s v - C%)(1 + 7s)]/2, the matrix 
elements of the operators are given as 



(<pK v \Of 



0\r^Z-FCNC\ 



Z-FCNC\ 
2 



B°) 



B°) = 



+ 



3 2 

4 (Cy + C%) 

3 2 



X 



nc v - 

3 2 



f7 s 



The values for C v and are taken as 



Cy 



Thus the amplitude for B 

as 



1 2 ■ 2n 

~2 + 3 Sm 



= -- 



(35) 



(36) 



<f)K arising from the Z mediated FCNC tree diagram is given 
G F 4 



<I>K, 



) = -^-(-l + sm 2 e w )U sb X . 



(37) 



Using the experimental upper limit Br(B -> X s l + l-) < 4.2 x KT 5 [36], in Ref. [35] 
the bound on \Uf, s \ is found to be \Ub s \ < 2 x 10~ 3 . Recently Belle Collaboration [37] has 
measured the branching ratio for the process B — > X s l + l~ as 



Br(B -> X,/+/-) = (6.1 ± 1.4li;f) x 10" 6 
Using the above result one can obtain the value [35, 38] 



0.06 ±0.03 , 



(38) 
(39) 



\Y (x t ) \T + C U2Z U bs \ 
where all the parameters in (39) are given in [35]. Thus one obtains the value of Ub s as 

\U bs \ ~ 1 x 10~ 3 (40) 



Now using sin 2 ^=0.23, we find 



r NP ~ 0.58 



(41) 



The variation of S^k with respect to strong phase 5np and weak phase 9^p according to 
Eq. (16) in VLDQ model is shown in Figure-2 and the variation of branching ratio (15) is 
shown in Fig-3. It can be seen from the figures that the observed asymmetry S^k and the 
branching ratio can be easily accommodated in this model. 



4 CP Violation in B — > rfKs process 

At this stage we are in a position to test, as mentioned earlier, whether the above two 
models (Model-Ill of THDM and VLDQ) can accommodate the result for another similar 
mode, which seems to be in agreement with the SM. In doing so, now we consider the 
B° — > rj'K process. 
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0np 90 

18 0° 



Figure 2: 3DPlot of S^k versus the weak phase 6^p and strong phase 5np (in degrees) for 
r NP = 0.58 

4.1 Contributions from SM and THDM 

In the SM, in addition to b —>■ sqq ( q = (u,d,s)) penguins, the B° — ► rj'Kg process also 
receives a small contribution from color suppressed b — > wws tree diagram. We first find out 
the standard model contribution. The matrix element in the SM is given as 



A(B^ V 'K°) = ^ 



10 



i=l 



- V tb VZJ2Ct"(rj'K \O i \B°) 



i=3 



(42) 



where Oi t 2 are the tree and 03_6(7-io) are the QCD (electroweak) penguin operators. The 
matrix elements of these operators are given in the factorization approximation as [14] 



{•n'K \O x \B ) = -X 2 ( V 'K \O 2 \B ) = X 2 

o 

( V 'K°\O 3 \B ) = ~Jf! + 2X 2 + 
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50 100 150 200 250 300 350 

0NP 



Figure 3: Branching ratio of B — > (f)K° process (in units of 1CT 5 ) versus the phase <p NP 
(in degree). The horizontal solid line is the central experimental value whereas the dashed 
horizontal lines denote the error limits. 



where 



X 1 
X 2 

x 3 
Ri 



(V'K° 
(rj'K 
(V'K° 
(rj'K 
(rj'K 



Ot\B°) 
5 \B°) 
6 \B°) 
7 \B°) 
8 \B°) 
9 \B°) 
O w \B°) 



2 4 
Xi -\- — Xo ~\~ —X 3 

3 3 



Ri 



X\ — 2X2 



R\X\ — -X 2 



X, 



3 

R2 ) X 3 



R\Xi 



y^o 

3 

R\Xi — ^ 

3 3 3 

3 3 



X„i 





R 2 \ 


x 3 








3 j 




(5- 


R 2/ 


)x 3 









= I 



m 2 B -m 2 v ,)F B ^(m 2 K0 )-^f K 
i(m 2 B -m 2 KO )F B - K 



2 \ x v > , 



i(m 2 B - m 2 K0 )Fr K {rn 2 )Y v ,J2p K - ft 



2m 2 K0 



(m b - m d )(m s + m d ) 



R 2 



2 (2m 



K° 



mz 



(m b - m s )(m s + m s ) 



(43) 



(44) 
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X v / = 0.57 and Y v i = 0.82 are the mixing parameters of the uu + dd and ss components in 
the rj' meson [39], which correspond to 9 P = —20°. Thus the amplitude is given as 



A{B Q ^r 1 'K°) = 

2 



+ (^2(a 3 - as) - i(a 7 - a 9 )^ X 2 

+ (a 3 + a 4 - a 5 + ^(a 7 - a 9 - ai ) + (a 6 - #2) ^3} 



(45) 



The decay width can be given by 

T(B° -> 7/K ) = -M-|A(£° -> V^°)| 2 (46) 
87rm^ 

Using Fo B_>7r ^(m|-o) = 0.335, fx^) — 0.16(0.13) GeV, the quark masses as (m d , m s , m b ) = 
(0.0076, 0.122, 4.88) GeV and the values of the coefficients a^s for N c = 3 from Ref. [27] 
we obtain the branching ratio in the standard model as 

BR(B° -> r]'K°)\sM = 3.24 x 10" 5 . (47) 

which is slightly less than the current experimental data [28] 

BR(B° -> r]'K°) = (5.8ll;|) x 10~ 5 . (48) 

Now we consider the contributions arising from THDM. As discussed earlier in this case 
due to the presence of new charged Higgs penguin diagrams, the values of the effective Wilson 
coefficients a^s get modified. Again substituting their values from [27] in Eq. (45), we obtain 
the transition amplitude as 

A THDM( B _^ ^0) = A SM^ + Q 27 e i(e NP+ 8 NP )^ _ (4g) 

Now taking r NP = 0.27 and varying the weak phase Onp = {—^, ^} and strong phase 5^p = 
{0, 2n} we can see from Figure-4, that the observed value of S^k can be accommodated in 
the THDM. Furthermore, the observed branching ratio can also be explained in this model as 
seen from from Figure-5. If we take a crude assumption that the THDM and SM amplitudes 
interfere constructively, the maximum value of branching ratio is found to be 

BR THDM (B° -> rj'K ) = 5.22 x 10" 5 , (50) 

which lies within the present experimental limits [28]. 
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Figure 4: 3DPlot of S v 'k versus the weak phase O^p and strong phase 5np (in degrees) for 
r*ro = (127 



4.2 Contributions from VLDQ Model 

Now we consider the contributions arising from extra vector like down quark model. In this 
case the B° — > i] f K° process proceeds through both color allowed and color suppressed tree 
level Z-mediated FCNC diagrams. The corresponding operators are given as 



Of 



Z-FCNC 



(^±^) (1 _ 75) + (^ 



n Z-FCNC 
w 2 



(Cy + 



75) + 



(Cy - C%) 



(1 + 75) \qp] (51) 



Using Fierz transformation and equation of motion the matrix elements of these operators 
are given as 

( v 'k \o^- fcnc \b ) = (C^ + C d A )X 2 + C A X 3 



+ 



6 



(C v {l+R 1 ) + C d A {l-R 1 ))x 1 
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Figure 5: Branching ratio of B° — > r/'K (in units of 1CT 5 ) process in THDM versus the 
phase 4>np (in degree). The horizontal solid line is the central experimental value whereas 
the dashed horizontal lines denote the error limits. 



( V 'K»\0 : 



+ (c v (l + R 2 ) + C A (l-R 2 )^X 3 



l -(C u A + C d A )X 2 + l -C s A X, 



+ 



(c d (l + R 1 ) + C d A (l-R 1 ))x 1 



+ (C V (1 + R 2 ) + C S A (1-R 2 ))X : 



So the amplitude for B° — > rj'K in VLDQ model is given as 



A VLD Q(B°^r)'K ) = ^U sb ^{(C u A + C d A )X 2 + C A X 3 } 

+ ^( y C d (l + R l ) + C d A (l-R 1 ))x 1 
+ (c v (l + R 2 ) + C s A (l-R 2 f)X 3 } 



Substituting the values of C v , A -> as 

1 



a 



V 



C 



(s,d) 
V 



we find 



- sin 2 6 W , 

1 2 • 2n 

- + -sin 9 W , 



r NP ~ 0.72 



1 

2 ' 



1 

2 



(52) 



(53) 



(54) 



(55) 
(56) 
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The variation of SL'ifO and the branching ratio according to Eqs. (16) and (15)in the vector 
like down quark model are shown in Fig-6 and Fig-7. It can be seen that the observed 
asymmetry and branching ratio for B° — > rj'K mode can be easily accommodated in this 
model. 




Figure 6: 3DPlot of S^kq versus the weak phase 9np and strong phase 5np (in degrees) for 
r NP = 0.72 



5 Conclusions 

To summarize, the time dependent CP asymmetry measurements in B —>■ <pK s gives sin 2/3, 
which is 2.7a deviation from the corresponding value in B —>■ ipKg. According to the SM 
expectation they should measure the same. Unlike the B — > if)Ks, which is a tree level 
process, B — > <pKs occurs at one loop level, which allows room for new physics to play an 
important role. In this paper, we have explored two simple beyond SM scenarios, the two 
Higgs doublet model (model III) and model with extra vector like down quark. We found 
that model III of THDM is unable to explain, whereas vector like down quark model can 
easily explain the result. 
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50 100 150 200 250 300 350 

0NP 

Figure 7: Branching ratio of B° — > rj'K process (in units of 10~ 5 ) in VLDQ model versus the 
phase 0ivp ( m degree). The horizontal solid line is the central experimental value whereas 
the dashed horizontal lines denote the error limits. 

It is important to note here that any new physics scenario that explains the <j)K s discrep- 
ancy must also explain another similar two-body decay B — > i]'K s , which is also expected to 
give the same value as of ipK s or <pK s , i.e., sin 2(3. In doing so it will be easy to rule out or 
narrow down the various NP scenarios. We found that both the models (THDM-model-III 
and VLDQ) can explain the r\'K s result. This in turn gives us the clue that the VLDQ 
model may possibly be a strong contender for the NP effects responsible in B — > <pK s . It 
is worthwhile to emphasize that various supersymmetric models (as can be found in the 
literature) can explain the (f>K s discrepancy. But apart from [13, 14] none of the scenarios 
so far explained the simultaneous explanation of 4>K S and r]'K s . On the other hand, our 
findings indicate that the simple non-supersymmetric extension of the SM in terms of the 
matter content should not be ignored for possible NP candidature. Regardless of the sources 
of NP, if in future the 4>K S result continues to be different from the SM expectation, then it 
will certainly establish the presence of NP. 
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